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Abstract. Let u{t) = —Fx{t) be the optimal control of the open-loop system x'{t) = Ax{t) + 
Bu{t) in a linear quadratic optimization problem. By using different complex variable argu- 
ments, we give several lower and upper estimates of the exponential decay rate of the closed-loop 
system a;' (t) = [A — BF)x{t). Main attention is given to the case of a skew-Hermitian matrix A. 
Given an operator A, for a class of cases, we find a matrix B that provides an almost optimal 
decay rate. 

We also show how our results can be applied to the problem of optimizing the decay rate 
for a large finite collection of control systems {A,B-j), j — 1,. . . ,N, and illustrate this on an 
example of a concrete mechanical system. Some open questions are listed. 



1. Introduction 

It is well-known that in many practical problems, an engineer has to optimize, in one or 
another sense, several performance parameters of a control system. The Linear Quadratic Opti- 
mal Regulator (LQR) problem searches a stabilizing feedback which optimizes some associated 
quadratic cost functional. Another important characteristic of stabilization is the exponential 
decay rate of the resulting closed-loop system. The main question we address in this article is 
to study in which situations the LQR provides good decay rates of the closed-loop system. 

Recall that the standard Linear Quadratic Optimal Regulator problem concerns the dynamic 
system of the form 

(1) x'{t) = Ax{t) + Bu{t), x{0) = xo. 
The problem is to minimize the cost functional 

/•oo 

(2) r{xo)= I x{tyQx{t) + u(t)*Ru{t)dt. 

Jo 

Here x{t) G is the state of the system and u G Lf^^(^[0, +qo),C"^) is a control function. 
Matrices A, B, R, Q are complex and have suitable sizes. We assume that R and Q are positive 
definite. We are specially interested in the case when the dimension m of the control u{t) is less 
than n, the dimension of the state x{t). 

As is well-known (see [19], [33]), the solution to the LQR problem is unique and the function 
u{t), for which the minimum of the cost functional is attained is given by the feedback function 
u{t) = —Fx{t) = —Fe^'^~^^^^XQ, where F = R~^B*X is the feedback matrix and X is any 
nonnegative solution of the continuous Algebraic Riccati Equation 

(3) XBR-^B*X -XA- A*X - Q = 0. 

This solution X is unique and positive definite, and the minimum cost functional is given by 
J(xo) = XqXxq- It is also notable that the feedback matrix F does not depend on xq. The 
closed-loop system is 

x'{t) = AcUoopX{t), 
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def 

where j4ci.ioop = A — BF is stable, that is, its spectrum crd.ioop lies in the open left half-plane 
C_. We denote by || • || the euclidean norm of vectors in C'^ and the induced norm of matrices. 

The linear quadratic problem is one of the most widespread methods for stabilizing systems. 
In this work, we give various estimates of the quality of this stabilization in terms of the geometry 
of the spectrum of the open- loop system matrix A and the characteristics of B. We remark that 
the pole placement problem is known to be very ill conditioned for control systems of large size 
and that the linear quadratic stabilization is one of the methods for overcoming this difficulty. 
We refer to [I !] for theoretical results and for a discussion (see especially Section 4). 

The exponential decay rate of the closed-loop system is given by 

(4) 7decay B) = TOITV [\R(iv\: U £ fJcLloop} • 

It is well-known that 

7decay = sup{e > : yxo3K = K{e,xo) : \\x{t)\\ < Ke-'\yt > 0}. 

Hence 7decay can be seen as a characteristic of the quality of the LQ control for large times t. 
The LQ regulator can be considered to be good in this sense if 7decay is big. 

The main results of this article concern upper and lower estimates of 7decay This is done under 
the assumption that the matrix A is skew-Hermitian: A* = —A (that is, iA is Hermitian). This 
assumption just means that under the absence of control {u{t) = 0), the energy ||x(t)|p is 
conserved. Notice that if an open-loop linear system models a mechanical (or electrical) system 
where the energy is conserved, then we are in this situation. 

We also will assume that 

Q = I, R = I. 

The assumption about Q is rather natural in view of the above remark on the conservation 
of energy, and the case of general R > reduces to the case when R = I hy making a linear 
substitution u{t) = R^/'^u{t) in (1). 

As we show, the upper and lower estimates of 7decay we give permit one to compare the 
performance of the LQ optimal regulators of control systems [A,Bj), in which A is fixed and 
there are several possibilities for the matrix B. 

We are not aware of any previous work estimating 7decay for LQ optimal regulators. Other 
measures of the quality of control have been studied already. Among the most popular of them 
are the eigenvalues of X, \\X\\, trace X and detX. Since J{xq) = XqXxq, these measures are 
tightly related to the cost of the stabilized system. 

Indeed, ||X|| has the sense of the worst case performance of the cost functional, for xq of fixed 
norm: 

||X|| = max J(xo). 
ll^o||=i 

Similarly, trace X is the average value of J{xq) when xq ranges over the unit sphere. The 
larger is any of these measures of quality of the control, the worse is the LQ stabilization. 

Estimates for all these measures are well know. See for instance the reviews by Mori and 
Derese [27], and Kwon, Moon and Ahn [IS], the papers [17], [26], [34] and recent papers [ ], [ ], 
[20], [21], [22]. 

We observe that a lower estimate of 7decay could be deduced from an upper estimate for ||X|| 
and the following easy relationship: 

1 

(5) 7decay > ^-jj^ • 

This last inequality is true because for any v G (Tci.ioop; if [A — BF)xo = uxq and ||xo|| = 1, then 

/■oo ^ 

\\X\\ > {Xxo,xo) > / \\x{t)fdt 







2 Rez/ 



Several works give upper bounds for \\X\\, but we are not aware of any estimate that would 
apply to the case when m < n and A is skew-Hermitian. The upper bounds given in these works 
either require that A + A* < or that BB* is invertible. 
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Notice that (5) shows that whenever the stabihzation is bad in terms of the parameter 7decay) 
||X|| also is large. 
We put 

a{A) = {iXi, . . .,iXn} 
(where Xj £ M) and assume throughout the whole article that 

(6) Ai < A2 < • • • < A„. 

Our estimates depend on the following numbers. The characteristic 

(7) 6{A) = min \Xj - Xk\ 

gives the minimal separation of eigenvalues. We will write just 6 when the dependence on A is 
clear enough. For a fixed index k, we put 

(8) 6k = min \Xj - A^l, 

which denotes the separation of the eigenvalue iX^ of A from the rest. The number 

(9) A(^) = max \Xj - Xk\ = A„ - Ai 

will also be used. 

The skew-Hermitian matrix A can be diagonalized: 

(10) Avj = iXjVj, 
where {vj} (1 < j < n) is an orthonormal basis of C". Put 

(11) bj = B*Vj. 
One of our main results can be stated as follows. 



Theorem 1. Put 

WhW 

l<fc<n^(l + 2 



Then the following statements hold. 

(1) The eigenvalues Vj of the closed-loop system lie in the box [— ||i?||, — ^est) x [Ai,A„]. 

(2) // moreover, m < n and the smallest singular value am of B satisfies am > 2\/2A, 
then exactly m eigenvalues lie in box [— — ^cJm] x [Ai,A„], and the other n — m 
eigenvalues lie in the box, (— \/3A, — £est) x [Ai,A„]. 

(3) In the case m = 1, the bound i^st in the above assertions can be improved by substituting 
it by a larger number 

(13) ^- = ,min 

In particular, it follows from this theorem that 

(14) Tdecay > 4st (idecay > ^Ist ^r 171 = 1). 

If A has multiple eigenvalues, we put £est = ^Ist = 0. It follows from the proof of this theorem 
that all its statements remain true in this case. 

Theorem 2 below gives a more detailed information about the location of the spectrum of the 
closed-loop system. 

Notice that the appearance of the norms of vectors bk = B*Vk in this estimate is very natural. 
In fact, the quantity 

do(A^) = min||6fe|| 

k 

can be taken for a kind of measure of controllability of the system x = Ax + Bu. In the case 
when all eigenvalues of A are distinct, the system is controllable if and only if min^ \\bk\\ > 0. 

We remark that iim = n and for some fixed (3 > one can freely choose B with \\B\\ = /3, then 
an optimal control with the best possible 7decay can be given easily. If BB* = /3^/ (for example 
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take B = pi), then the solution to the associated continuous Algebraic Riccati Equation is 
X = I3^^I. Hence, the closed-loop system matrix is ^ci.ioop = A — f3I, and one can readily 
compute its eigenvalues. If follows that in this case, in the bound 7decay < 11-^11) which follows 
from Theorem 1, the equality is attained. 

We also observe that the case m > n can be reduced to m < n. In fact, the optimal feedback 
u{t) = —B*Xx{t) ranges over the space Rani?*. Therefore the linear quadratic problem for the 
pair {A, B) reduces to the same problem for the pair (^4, B\ Ran S*) ; notice that dim Ran B* < n. 
After this reduction, in place of B, we get the operator B\ Rani?*, which has trivial kernel. 

For this reason, we will assume throughout the paper that 

m < n and ker B = 0. 

Let us briefly overview the contents of the article by sections. Section 2 is devoted to the 
proof of Theorem 2, which implies Theorem 1 above. 

In Section 3, we show that if the minimal separation d{A), defined in (7), is rather big in 
comparison with ||i?||, then the closed- loop eigenvalues of the system can be located with good 
precision, which gives nice two-sided estimates of 7decay In particular. Corollary 11 shows that 
if ||i?||/(5(^) is rather small, then 7decay is comparable with min \\bk\\- In many problems of the 
design of optimal controllers, the matrix B can be changed, up so some extent. In this section, 
for a given A, we find a "suboptimal" matrix B among all matrices with a fixed norm, which is 
supposed to be small. (See Theorem 12 and Corollary 13.) 

For l<A;<n — Iwe define 

(15) Afc =^ min |Aj+fc-Aj|. 

l<i<n— fe 

Observe that with this notation, 5 and A defined in equations (7) and (9) are 5 = Ai and 
A = A„_i. 

In Section 4, Theorem 14, we give an estimate of 7decay in terms of A^ (recall that m is 
the dimension of u{t)). For m > 1, this estimate may be much better than the estimate of 
Theorem 1 if some of eigenvalues of A are close to each other or coincide. 

Section 5 contains a brief account of all our estimates of 7decay In Section 6, some numerical 
examples that illustrate these estimates are given. In Subsection 6.1, we give an example in low 
dimension, which illustrates how our estimates compare in different cases. In Subsection 6.2, 
we discuss the problem of optimizing 7decay among a finite family {A,Bj), with a fixed system 
matrix A and different possible choices for the control matrix B. We give an algorithm which 
uses our estimates to reduce the number of computations needed in the search. We illustrate 
this algorithm with a simple mechanical system. 

In Section 7, we list some open questions, and in Section 8, we list the conclusions of this 
article. 

In what follows, we use the notation for the Frobenius norm of a matrix Y. It is given 

by 

= trace(y*y); 

this formula applies to rectangular matrices as well. 

It is worth noticing that for general pairs of matrices {A,B), Eising in [ ] introduced the 
so-called "distance to uncontrollability" , given by 

duc{A,B) = mf{ (\\6Afp + \\SBfpy^ : SA G C"^", 6B e C"^"^, 

(16) ^ . 

{A + 6A,B + 6B) uncontrollable 

He proved that 

duc(.A,B) = mm a^in{[A- XI, B]), 

AG*!- 

where <Tmin stands for the minimal singular value. Denote by rfc(cj(A)) the radius of the smallest 
disk containing at least k points of cr{A). It is not difficult to show that for any normal matrix 
A and for any B such that m < n, one has an estimate 

duc{A, B) < min {do{A,B),r^+i{a{A))). 
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Figure 1. Location of the eigenvalues of ^ci.ioop- 

If iA is Hermitian, then rm+i{o'{A)) = ^A^- Estimates for the quantity duc{A, B) and methods 
for its computation have been studied further in numerous works, see [2], [8], [10], [13], [16], 
[31] and references therein. Related characteristics were studied in the works [32], [15] and 
others. Numerical methods for solving matrix algebraic Riccati equations, in fact, have been 
much investigated; see [4], [28] and books by Sima [30] and Datta [7]. We refer to [20] and to 
[3] for some other interesting aspects of the linear quadratic problem. 



2. The main result on the location of closed-loop eigenvalues 
The spectral theorem yields the decomposition 

n 

(17) A = Y,i^3V3y*v 

i=i 

where the eigenvalues i\j of A are assumed to satisfy (6) and the eigenvectors Vj form an 
orthonormal basis of C" (see (10)). Moreover, B and B* decompose as 



(18) B = Y,v,b] B* = ^b, 

i=i i=i 

where the bj^s have been defined in (11). 

For any index k, 1 < k < n — 1, we consider the closed right triangle in C with vertices at 
the points 

iXk, ^Xk+l, ^ H^ . 

All these triangles lie in the half-plane Rez < (see Figure 1). For any fc, 1 < A; < n, we put 

(19) Pk = z zm^ ' Pk ~ 



Our next goal is to prove the following result. 

Theorem 2. (1) The eigenvalues Vj of the closed-loop system lie in the box [— ||i?||,0) x 
[Ai, An], outside the triangles and outside the closed disks centered in iX^ of radii pk, 
given by (19). 
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(2) // moreover, m < n and the smallest singular value am of B satisfies am > 2^/2A, then 
exactly m eigenvalues of the closed-loop system lie in box 

[-\\B\\,-^am] X [Ai, An], 
and the other n — m eigenvalues lie in the box 

(-\/3A,0) X [Ai,A„]. 

(3) In the case m = 1, the assertion of (1) holds for disks with the same centra and larger 
radii p\, instead of pk- 



Notice that lest = niirifc /)fc/v2, £g^j = miiifc pyV'^ (see Figure 1). Therefore Theorem 1 is an 
immediate consequence of the above theorem. 

Before proving Theorem 2, we need some prehminaries and several lemmas. 

2.1. The function $ and its zeros. The rational matrix function defined as 

(20) ^{z) = I - B*{zl - Ay^B. 

is important in the control system theory. It is known that $ factorizes as 

<^{z) = M{-z)*M{z), 

where 

M{z) = B*X{zI - Ay'^B + I. 

The theory also shows that 

M(z)-^ = -B*X{zI -{A- BB*X))-^B + /. 

See, for instance, the book by Zhou, Doyle and Glover [ , chapter 13.4] for a proof of this 
factorization. 

Hence, the eigenvalues of A—BB*X are poles of M{z)'^ in the sense that if zq is an eigenvalue 
of ^ — BB*X then det(M(2;o)~^) = oo. It follows from the factorization of ^{z) that the zeros 
of ^{z) (in the sense that det ^{z) = 0) are 

{z G C : det $(z) = 0} = {2 G C : 2; G crd.ioop or - z G cJci.ioop}- 

Definition. Let ^{z) be as in (20) and G C such that det ^{zq) = 0. If Rezo < 0, then zq is 
called a stable zero of ^{z). If Rezo > 0, then zq is called an anti-stable zero of ^{z). 

So the stable zeros of $(z) are exactly the eigenvalues of ^ci.ioop* 

The function <I> will be very useful to make estimations of the cost characteristic 7decay The 
relation between 7decay and ^[z) is 

(21) 7decay(A,5) = min{|Rez| : det$(z) = 0}. 
If we define 

then, for fixed z G C, / is holomorphic in A on C \ {z} and hence f{A]z) is well defined if 
z ^ (y{A). We can write 

^z) = I-B*f{A-z)B. 

Using (17) and (18), we get 

n 

(23) 

An important remark is that ^{z) is Hermitian and positive along the imaginary axis where 
it is defined. Indeed, we have A = iAq, where ^0 = ^o- t G M, t ^ Xj, then 

$(it) = I- B*{itl - A)'^B = I + B*{tl- Aoy^B > I, 

because B*{tl — Aq)^'^B is Hermitian and positive. 
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Lemma 3. The zeros of ^ lie in the box in the complex plane given by |Re2;| < \\B\\, Ai < 
Imz < Xn- 

Proof. Recall that the real and imaginary parts of an operator T are defined by 

Rer = ^(r + r*) imr = ^(r-r*). 

Put 

h{X; z) = \{f{\- z) + /(-A; z)), g{\- z) = i(/(A; z) - /(-A; z)). 
Then h and g are meromorphic in A on the whole plane. It is easy to see that 
(24) Re^{z) = I - B*h{A;z)B, lm^{z) = -B* g{A; z)B. 

z = x + iy , \ €z M, a. direct computation shows that 

x'-{y- A)2 



(25) h{iX;z) 

(26) g{iX;z) 



(x2 + (y-A)2)2^ 

-My - A) 



(x2 + (y-A)2)2- 

First we show that if | Rezj > ||i?|| then Ke^{z) > so z is not a zero of ^(z). Let G 
with ll^ll = 1. Then 

{Re^{z)(„^) = 1- {h{A;z)B(„B^) >l-\\Bf max h{X;z). 

Aeo-(A) 

Now, using (25), if | Rez| > ||i?||, it follows 

1 



max/i(A; z) < ,, ^,,„ 
Ae« ||-B||2 

and therefore (Re<^(z)^,0 > for these z. 

Now observe that if either Imz < Ai or Imz > A^ then (26) shows that g{X; z) has constant 
sign for all A G <y{A) and therefore Im<I>(z) is either postive or negative (since we may assume 
keri? = 0) so that z is not a zero of ^{z). □ 

In Lemma 3 we have seen that the zeros of ^{z) cannot be too far from the imaginary axis. 
The next two lemmas imply that the zeros cannot be too close to the imaginary axis. 

Lemma 4. Define the angles 

2lfc = {z G C_ : ^ < arg(2; - iXk) <-^]-, 1 < k < n. 
If z is in the left half-plane, but does not belong to the union of these angles, then Re<I>(z) > I. 
Proof. Put z = X + iy. It follows from the hypothesis on z that 

x^-{y-Xjf<0, j = l,...,n. 

Defining h(X; z) as in the proof of Lemma 3 and using (25), we get that h{A; z) < 0, so that 
Re^>(z)>/. □ 

It follows from the above two lemmas that the stable zeros of $ lie in the band Ai < Im z < Xn 
and outside the triangles Ti, . . . ,T„_i. 

Lemma 5. ^ has no zeros in the disks D(iXk, Pk), I < k < n. 

Proof. If pfc = then the lemma is vacuously true for the corresponding k. Hence, assume 
Pk > 0. Suppose Rez < and \z — iXk\ < Pk for some k. Fix this index k. Observe that ^ is 
the length of the legs of one of the triangles T^, whose vertex is in iXk- Since 



yk , _ II D||2 X 2 ^ nil DII2 — 



w ^ \\hk\m ^ K 



;i + 2^^2 2II5II2 - 2' 

it follows that z belongs to D(iXk, 



V2' 
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It can be shown geometrically that the intersection of D{iXk, with D{iXj, (j ^ k; 
notice that we take the same radii) and with the left half-plane is either empty or is contained 
in one of the triangles T^. Therefore, if z G D{iXj, ^) for some j 7^ k, then z is inside of one of 
the triangles, and it has already been shown that then z will not be a zero of $. 

So let us assume that 

(27) \z — iAyfcl < , but \z — iXj\ > for j ^ k. 

V 2 v 2 

Put 

(see (17)). It follows from the first inequality in (27) that z ^ 2tj for j 7^ k. Hence Lemma 4, 
applied to the configuration of n — 1 points {iXi, 1X2, • • • , «A„} \ {iXk} on the imaginary axis, 
gives that ReC(z) > / (recall we have assumed keri? = 0). 
By (27), we also have 

||C(z)||<M^=l^l + ^. 

Next, let us show that the above properties ||C(2;)|| < M and ReC(z) > /, imply that 
ReC(z)-i > Indeed, take any ^ G with ||^|| = 1 and set ry = C-^{z)i. Then 

1 = ll^ll < ||C(z)||||r/|| so that ||r/|| > }j. Hence, 

Re{C{z)-^i , = Re(C(z)7?, r?) > ||r?f > ^ , 

and the inequality ReC(z)^^ > follows. 

Now suppose that ^{z)^ = for some fixed ^ G C™, ^ 7^ 0. Then, 

Since C{z) is invertible, we have 7^ 0. Multiply the above equality by (6|.^)^^6|.C(z)^^ and 
regroup terms to yield 

hlC{z)~^hk = {z-iXkf. 

Since ReC(z)~^ > we get 

|.-zA,|2 = |6^C(.)"i6,|>M!=^2, 

Therefore \z — iXk\ > p^, a contradiction. □ 
In the case m = 1, the above lemma can be strengthened. 



Lemma 6. If m = 1, then $ has no zeros in the disks D{iXk, p\). 



Proof. Assume that m = 1, ^[z) = and for some index k, z £ D{iXk,pl.)- Proceed as in the 
previous lemma to deduce that I2: — Aj| > ^ for j 7^ k. Now, since ^{z) = 0, we have 

{z - ix,Y j^^ {z - a,)2 

(notice that now hj are complex numbers). It follows that 



< 



so that \z — iXk\ > pj., a contradiction. □ 
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Lemma 7. Let am be the minimum singular value of B. If < ^Um, then exactly m of 
the stable zeros o/$(z) lie in the box given by 

— ||i?|| < Rez < — ^'^"i' Ai < Ima; < A„, 
and the n — m remaining stable zeros lie all in the box given by 

-VSA < Re z < 0, Ai < Im z < A^. 
In particular, no stable zero lies in the band Rez G {—^o'm, — \/3A]. 

Proof. The restriction to Ai < Imz < An comes from Lemma 3. To prove the statement about 
boxes, suppose that am satisfies the hypothesis given. 

Let Tc : [a, /3] — )• C be the closed positively oriented contour, traversing the boundary of the 
box [— c, c] X [— (i, d]. Since <I>(oo) = I, d can be chosen large enough so that all the eigenvalues 
of ^{z) are arbitrarily close to 1 when z is on the horizontal segments of Tc- We assume that 
Tc{a)=Tc{l5) = -c-id. 

Let 7c be the right vertical segment of Lc, going from c — id to c + id. We subdivide 7c into 
three segments, 

[c-id,c + i\i], [c + i\i,c + i\n], [c + i\n,c + id]. 
We will use expressions (24) for the real and imaginary parts of ^{z). First observe that if 
z E [c — id,c + zAi], then Im$(2;) < 0. Indeed, for these z, Rez > and Imz < Aj for all j. 
It follows that g{iXj;z) > and therefore lm.^{z) < (see (26)). Hence, all the eigenvalues of 
^(z) lie in the open lower half-plane. 

Similarly, if z G [c + iXn, c + id], one has Im$(z) > 0. Hence for these z, all the eigenvalues 
of ^{z) lie in the upper half- plane. 

Now we will show that if 

/fi 

(28) VSA < c < ^am 

then Re $(2;) < for z G [c + iXi, c + iA„]. Write z = c + iy, y £ [Xi, A„]. Then, using (25) and 

(28) , we get that for all j, 

hC\ c^-{y-X,f c^-A^ c^-c^/'i _ 3 J_ 

^(*Aj; Z) - ^^y_ > (^2 + A2)2 - (c2 + c2/3)2 ~ 8c2 ^ al' 

If ^ e C™ with ll^ll = 1, then 
(Re <^{z)^, = 1- {h{A; z)B^, B^) < 1 - { min h{X; z)) ■ < 1 - -^mW < 0, 

X<^a(A} a^ 

because \\B^\\ > amUW- Hence, Re<i>(z) < 0. 

Since ^{—z) = $(z)*, ^(z) behaves similarly on the left vertical segment of Fc. 

Now choose c satisfying (28) and study the winding number of det <I* o Fc around 0. The m 
eigenvalues of $(z), fi{z), . . . ,ipm{z), can be numbered so that ipj o Tc{t) are all continuous 
functions of the parameter t,t £ [a, /?]. Since det $ o Fc = {^fi o Fc) • {ip2 ° Fc) • . . . • {ipm ° Fc), it 
follows that 

m 

(29) index(det <I> o Fc) = ^ index((/7j o Fc). 

Let us calculate the winding number of the curves ipj o Fc. 

When z = Tc{t) is in the lower horizontal segment, ipj are all close to 1. Then, as z travels 
through 7c, (fj first are all in the lower half-plane, then go to the left half-plane and then to 
the upper half-plane. When z is in the upper horizontal segment, all the numbers (pj are again 
close to 1. It follows that by choosing d sufficiently large, we can make the winding number 
of each of the functions ipj to be arbitrarily close to —1 on each of two vertical segments of Fc 
and arbitrarily close to on the two horizontal parts of Fc. Since det <^ o Fc : [a, /3] — )• C is a 
closed curve, its winding number around is an integer. By (29), it is equal to —2m. Using the 
argument principle and the fact that det ^{z) has 2n poles counting multiplicities inside Fc, one 
gets that det$(z) has 2n — 2m zeros inside Fc. Hence, ^{z) has n — m stable zeros inside Fc- 



10 DANIEL ESTEVEZ AND DMITRY V. YAKUBOVICH 

Setting c = V^A, we see that there are n — m stable zeros inside (— -y/SA, 0) x [Ai, A„]. Letting 

c — )■ ^(Jm, we obtain that again $ has n — m stable zeros inside the box {—^am, 0) x [Ai, An]. 
The remaining m stable zeros must lie all outside this box, and by Lemma 3, they belong to the 
boxto X [Ai,A„]. □ 

Proof of Theorem 2. All the statements of this theorem follow from Lemmas 3-7. □ 

Proof of Theorem 1. As we already pointed out just after the statement of Theorem 2, Theorem 
1 is its direct consequence. □ 

Using Theorem 2, we can provide upper bounds for the value of 7decay 

Corollary 8. The following upper boufid (ilways holds fov the value of^f^Qc^y^ 

(30) 7decay(^,5) < 

If in addition, am, the smallest singular value of B, satisfies am > 2^/2 A and m <n, then 

(31) 7decay(^,5) < \^A. 

Proof. The first bound comes from lemma 3. Under the conditions of the second bound, using 
lemma 7 it follows that ^{z) has at least one zero on (— \/3A, 0) x [Ai, A„]. □ 

We remark that, as we already mentioned before, Theorem 1 also implies lower estimates for 

7decay i 

namely 7decay > ^est and 7dccay > ^Lt for = 1. 
The following upper bound holds for l^st- 

, . llbA.11 5^ 
lest = mm — ^TTTT^ < 

Indeed, the inequality 



i<fc<'t ^(1 + 2^^^'! 4 \ n\\B\ 



ii&^ii<V^i^i 



(see the proof of Theorem 12) implies that we get 

\/2 . \\bk\\5l ^ V2 6' 



iest < — mm ^ < 

4 k \\B\\^ 4 y n \\B\ 

Similarly, i^^f < 5/{2y/n). One can guess that a matrix B in which all are as big as 
possible can be used to ensure a nearly optimal stabilization of the system. A matrix with these 
characteristics will be given below in Theorem 12. 

3. The estimate of 7decay IN THE CASE OF A SUFFICIENT SEPARATION OF THE OPEN-LOOP 

SPECTRUM 

Here we will assume that the minimal separation 5{A), defined in (7), is rather big in com- 
parison with ||i?||. We will use the following analogue of Rouche's theorem for matrix- valued 
functions. 

Lemma 9. Let F[z) and G{z) be meromorphic functions on some open subset D C C, whose 
values are m x m complex matrices. Let ^ be a closed curve in D such that det F{z) has no 
poles or zeros on 7. // ||F~^(z)G(2)|| < 1 for all 2; G 7, then the scalar functions detF(z) and 
det(F(z) + G{z)) have the same winding number along 7. 

This lemma is known. See, for instance, [12, Theorem 2.2] (for operator-valued functions) or 
[25]. 

The next theorem locates the points of the closed-loop spectrum inside disks of radii such 
that Tfc — 7- as " ^j' — ^ 0. Recall that the zeros of <I> which lie in the left half-plane coincide 
with the eigenvalues of the closed-loop system (see Subsection 2.1). 

Theorem 10. Suppose k is an index such that < . Then <I> has exactly one zero 

^ l|2 ,, ,, 



in the open disk of centre Zk = — \\bk\\ + i^k and radius r^ — ^2 



2 
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Proof. Observe that < and consider the contour 

7 = {zfc + rfcC*'' : < 6* < 27r} 

and the functions 



F{z) 



G{z) = ^{z)-F{z). 

These functions are holomorphic on 7 and its interior. We will prove that ||G(2;)|| < ||~-^ 
for z G 7, so that we can use the above version of Rouche's theorem. 

First observe that F(z) is normal, so that ||~"^ = mmx^(j(^F(z)) \M- The spectrum of 

F(z) can be computed easily: 

Take any z such that \z — Zk\ = and put z = Zk + r^e*^. Notice that < implies 
We'' -2\\hk\\ \ > \ \\hk\\ - rue'^l We get 



1 



{z - iXk)' 



\bk\ 



{\\bk\\-rue''Y 
rk\rke'' - 2 
\\\hk\\-rke''\ 



2rfc||6fc||e^' 



> 



rk 



rke 



> 



ie\2 



rk 



> 



rk 



\hk\\+rk 2\\bk\ 



and so llFfz) 



-iii-i 



— 21161.11 ' -^^^^^ ^'^^^ suffice to prove that ||G(z)|| < 2ifftriT- (^'^) 



Then, observe that the condition ^^^L < implies < (\/2 — l) 5k since 

rk < \\bk\\ < \\B\\ < ^^=^<5fc = {V2-l)5k. 



Now we have for z G 7 and for all j ^ k 



\z — iXj\ 



Tke 



id 



\\bk\\ + iXk - iXj\ > \iXk - iXj - \\bk\ 



rk 



> \Xk - Xj\ -rk>6k-{V2- l)Sk = (2 - V2)6k. 



Hence it follows that 



\\Giz)\\ < 



\B\ 



rk 



2Y51 nbk\ 



< \\F{zy 



-iii-i 



So, by Rouche's theorem, det F{z) and det(-F(z) + G{z)) have the same number of zeros inside 
7. The only zeros of detF(z) are Zk and —Zk- Since 7 lies completely in the left half-plane, 
det F{z) has exactly one zero inside 7. Therefore det $(2;) = det(F(2;) + G{z)) has exactly one 
zero inside 7. □ 



dcf 2||i?| 



Corollary 11. Set 

(32) ipk"= — '^^^ ^, l<k<n. 

Suppose (fik < 1 for at least one index k. Put 

r+ = min{(l + (/JA..)||6fc|| : v'fc < 1}, T 
Then 7decay < r+. If moreover ipk < 1 for all k, then 

T- < Tdecay < r+. 



min(l - Lpk)\\bk\\- 
k 
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Proof. If (fk < 1 for some k, the preceding theorem shows that some eigenvalue of the closed- 
loop system satisfies — Keu < (1 + and the upper bound follows. If ipk < 1 for all k 
and v is any eigenvalue of the closed-loop system, then — Rei^ > (1 — </5fe)||?)fcl| foi' some k, so 
that the lower bound follows. □ 

Using Theorem 10, when S{A) is sufficiently large, we can give a matrix B, in a sense close 
to optimal. 

Theorem 12. Suppose m<n. Let w he the primitive n-th root of 1 given by 

w = e " . 

Let ^ > 0. 

Let the matrix B he represented in the orthonormal basis given by {vj}, the eigenvectors of 
A, as 



(33) 



B 



(3 



w 



w 



0- 

1- 



w 



w 



0- 1 

1- 1 



w 



w 



0- (m-l) 

1- (m-l) 



yjin-iyo iijin-iyi ... y^(n-l)-(m-l) 

Then, \\B\\ = (3 and for any e > 0, there exists K > Q such that if 6 (A) > K, then 

( sup 7decay(-4,5)) - 7decay -B) < £■ 
\\B\\=^ 

Proof. First observe that B is related to the unitary Discrete Fourier Transform. If U £ c^x"- 
is the matrix of the unitary DFT, then 



B = (3U 



Ir, 







n—mxm 



It follows that ||i?|| = f3. Define bj in the same way as in as in (11), that is, put bj = B*Vj. 
Then \\bj \\ = /3y^ for all j = 1, ... , n. 

Let B be arbitrary with \\B\\ = /3. Let e > be given. Define K > from 

2 „o 4 (3^' 



Suppose (5^ > K^. Then ^ < "-^ , and the hypothesis of Theorem 10 is satisfied for any 
index k. We obtain disks of radii such that the zeros of ^ lie inside this disks. Now, 



max 



0^ 

2 - ^/2)2 (2 - V2y £ 



\\bk\\ < 



< 



13^ ^ 8 



{2 - 51''""'' " (2-^2)2^2 - (2-V2)2i^2 - 2' 
Since there is a zero of ^ in each of these disks of centre = — H^fcH + i^k-, we have 



min ||6fc|| 



< 7decay(^, B) < min + -• 



Notice that 



^\\bjf = \\B\\l <m\\B\\^ 



Hence, min 1 1 6,- 1 1 < B —. Therefore 



sup 7decay(^,^) < I^J h 

B|l=/3 V ^ 



Since B has 1 1 ft. 



/3 



for all j 



7decay 



m £ 



and the theorem follows. 



□ 
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Corollary 13. Assume that m < n. Let f3 > 0. If A is such that 

(34) S\A)> 

then the matrix B given in (33) satisfies 

(35) 7decay(^,^) > SUp 7decay(-4, 5). 

2 \\B\\=I5 

Proof. By applying Theorem 12 and its proof with e = one gets 

(2-^/2)2 V 

Therefore, if > then 

/rn e 2 ^ /rn 1 / ^ /rn e\ . 1 

\\B\ 

4. Estimates of decay rate in terms of for m > 1 
We begin with the foUowing remark. Let 9 > 0. Consider the following function 

f( )= ^"^ 

cT + a2 + ... + a'^+i' 

which is positive on (6*, +oo) and vanishes at 6 and at +oo. We denote by ao{6) the point in 
{6, +oo) where / takes its maximal value and by fi{6) = f[ao{6)) this maximal value. If we put 

m 

P{a) = a + a^ + --- + - ( ^(j + l)a^) ■ {a - 9) 

j=0 

then it is easy to see that ao{9) is the unique root of P in {9,+oo) (notice that P' < on 
i9,+oc)). 

In what follows, Pun{wi,...,wr} will stand for the orthogonal projection onto the linear span 
generated by vectors wi, . . . ,Wr- 

Theorem 14. Let 2 < m < n. Let vj be the eigenvectors of A (see (6) and (10)). Put 

Hfc = Pun{v,,...,v,+^.,} : ^ C", k = l,...,n-m + l. 
Assume that positive constants 7, K satisfy 

(i) \\B\\ < K; 

(ii) \\B*I[kd\\ > 7||nfcd|| for all d £ and for all k. 
Define Am from (15). Put ao = (To(^), 



7" 

(36) p=-mm|A„(2^-^-l) , 7^:^) 

Then all eigenvalues of the closed-loop system lie in the half-plane Rez < —p < 0. 

Notice that can be positive even in the case when some of the eigenvalues of A coincide. 
We do not exclude this case. 

The rest of this Section is devoted to the proof of this theorem. 

The plan of the proof is as follows. First we remark that it is easy to get from (i) and (ii) 
that K/7 > 1. Hence uo > 1. 

Fix some z = x + iy £ C_ such that —p < x < 0. We have to prove that z ^ Tcijoop- To do 
that, let us consider a reordering A^(i), At-(2)) • • • ; ^r(n) of the eigenvalues Ai, . . . , A„ of ^A such 
that 

(37) \Z - iXr{l)\ <\Z- iXr(2)\ < ■■■ <\Z- iXr{n)\- 
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Let US assume that 

(38) |z - A^(i)|^ < 

(if it is not true, then z ^ o"ci.ioop) due to Theorem 2). We will divide the spectrum (t(A) into 
two parts: 

fclose(^) = {^At-(i), . . . ,iA^(5)}, (Tfar(A) = {iA^(s_|.i), . . . ,iA^(„)}, 

where the index s will be elected according to Lemma 15 below. (Notice that the reordering 
(37) and this partitioning of cr(A) depend on the position of z.) Once this partition is chosen, 
we put 



Introduce the notation 



??close = \z — ^\{s)^ 1 ^far = 1^ " ^\{s+\)^ ■ 



(so that cj < x). 

We will say that adosei^) CLn-d cr/ar(^) are sufficiently separated (with respect to z) if 
(39) ryfar > 2x and ryclose < 



??far - 2x2 + X 

Inequality r/far > 2^^ implies that ?/far(w — x + 2^^ — r/far) < 0. By using the second inequality 
in (39), one gets that the sufficient separation implies the strict inequality r/ciose < ^?far- 
Before finishing the proof, we need three lemmas. 

Lemma 15. For any z such that —p < Rez < and (38) holds there exists an index s, 
1 < s < m, such that the corresponding parts afar{^) and adosei^) of the spectrum of A are 
sufficiently separated. 

Proof. Take some z that satisfies the hypotheses. Let A,-(i), . . . , A^(() be all point of the spectrum 
of A that satisfy 

r]j =^ |z — xAt-q^I^ < 2x'^. 

Since 

\x\< p< Am/2 

(see (36)), it follows that I < t < m and that r]t+i > 2x^ . Assume that the subdivision of the 
spectrum of A into two sufficiently separated parts is impossible. Then 

r\ ' — 2'3p' 

(40) ^^-^ 2:^:2 + x J- = 

We will prove that this leads to a contradiction. Put b = uj — x/o-q. Since uo > K'^ j^"^, we have 
< 5 < to! and ao = ^jz^. 
We prove that 

(41) m < ^ — 

do — i 

for J = t, . . . , m + 1 by induction in j. The induction base, j = t, follows from our assumptions. 
Assume that (41) holds for j = jo, t < jo < m. By using that |x| < /? and (36), we get 



2x^ <2p^ <^ p 



Therefore 

2 - 1 1 , X , 6 
Vjo < 2^' ^ <-{u--) = -<u;. 

(Tq — i Z (To Z 

We also have r?jo+i > 2^^. Hence (40) implies that 



2 Vio ^ 
Vjo+i - 2x < X- — — < r/j„- = r/joCJo. 



It follows that 
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This gives the induction step. Hence (41) holds for all j = t, . . . ,m + 1. In particular, 

r]m+i = \z- iK{m+l)\^ < 2x 



This gives a contradiction. Indeed, it follows that \r{i)^ • • ■ i \{m+i) contained in the interval 
[y- l,y + £], where I = {rj^+i - x2)V2, Then 

£2 < ^2 2 _Q 1 < 2 -2 1 < ^ 

~ V cJo-1 / V o-Q-l /~4 

(the last inequality is due to (36)). Hence < We get a contradiction to the definition of 

Am. □ 

Next, we take s = s{z) as in the above Lemma and put 

-Pclose = -fLm{?;^(i),...,D^(^)}) -Pfar = -fLin{D^(^+i) 

(recall that Vj is an eigenvector of A corresponding to iXj). Then 

(42) ^{z) = I - B*{z - A)-'^B = <i/ + T., 
where 

(43) * = / - B*Pi,,{z - A)-^Pi^,B, S = -B*P,iase{z - A)-^P,ioscB. 
Put 

W = B*Pt,„ 

then 

* = / - W{z - A)~^W* =1 + ^^, 

where 

^1 '^^^ -W{z- Ay^w*. 

Define a > from the equation 

(44) r?fer = (2 + a)x2. 
Then 

(45) (y - Afc)2 > (1 + a)x2 for all iXk G afar(^). 
Hence 

(46) {yl + fA)2 - x^lj RanPfar > ax^/| RanPfar. 
Lemma 16. Suppose \\B\\ < K. Then 

(47) Re^-^>/3-^I, 
w/iere 

13 = 1 + 



Proof. A calculation gives 

(48) Re^i = --W[{z- A)-^ + {z + A)~^]W* = WGW* 



where 



(y + iAf-x^I 

RanPfar. 



^^^^ ^ (z-^)2(z + ^)2 

We wish to prove that 

(50) + ^*-^ > 2/3-^1. 
First let us check the inequality 

(51) ^-^^ < ^(^ + ^'*) =/3Re^'. 
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Inequality (51) is obtained as follows: 

= + (7 + ^^) 
= / + 2Re^i + 

= / + 2WGW* + W{z + A)-'^W*W{z - A)-^W* 
<I + 2WGW* + K'^W{z + Ay^{z - Ay'^w* 

< / + (2 + — ^\WGW*. 

The last inequality is due to (46) and (49). By (48), this implies 

(52) </ + (! + /3)Re^'i = (1 + /3) Re ^ - /3/. 
Rewrite this inequality as 

or, equivalently, — ^^^J^-^H < This gives the inequality Re ^' < + ^~2^)^ = (^^- Then 
by (52), < 1(3 + l)Re^ - f3I < (3Re^, and we get (51). 

We also get that a{^) C {z : \z — ^^^| < ^y^}, which implies that ^ is invertible. So (50) 
follows immediately from (51). □ 

Lemma 17. Suppose 1 < s < m, is an m x m matrix satisfying Re^'^"'^ > fi^^I, V is an 
m X s matrix such that 

(53) \\Vd\\>^\\d\\, dea 

and A is an s X s invertible matrix. If P,^ are positive and 

(54) /3-V>l|A"^||, 
then the matrix ^ + yAl/* is invertible. 

Proof. Suppose it is not, 

(55) {"^ + VAV*)w = 
for some w € C", w 0. Put V*w = c G C*. One gets 

w = -^~^VAc, 
-V*'^-WAc = c. 

Hence 

(c,Ac) = -{V*'^-^VAc,Ac). 

Then, on one hand, one has 

\{c,Ac)\ = |(A-iAc,Ac)| < \\A-^\\ ■ \\Acf, 

and on the other, 

\{c,Ac)\ = \{V*^-^VAc,Ac)\ > Re(y*^-VAc,Ac) > l3-^\\VAcf > /J-^^IIAcf . 

These inequalities and (54) imply that c = V*w = 0. By (55), it follows that "^w = 0, so that 
w = 0, which contradicts to the choice of tt;. □ 

The end of the proof of Theorem 14- As before, we assume that some z = x + iy ^ C_ with 
—p < X < has been fixed. Lemma 15 gives us an index s, s < m, which defines a partition of 
a{A) into two sufficiently separated parts, (Tciose(^) and (Tfar(^). Define ^ and E from (43). 
By Lemma 16, (47) holds. Put 

V = B*\ Ran Pclosc , A = (z - A) | Ran Pdose 

and apply Lemma 17 to these two matrices and ^. Since r/ciose < Vfar, it follows that there is 
some k, l<k<n — m + 1, such that all the indices r(l), . . . ,r(s) are contained in the set 
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{k, k + 1, . . . , k+m — 1} (this is true even if A has multiple eigenvalues). Therefore, by hypothesis 
(ii) of the Theorem, V satisfies (53). By (44) and (39), one has 

n-l 2 2 '^•^^ 2 ^far ~ IIA-llI 
^ 7=7 2^^ =7 ^T^^ > ^close = A . 

Hence (54) holds. So <l>(2;) = + VAV* = ^ + is invertible, and therefore z ^ acUoop- This 
proves the Theorem. □ 

5. A BRIEF ACCOUNT OF OUR ESTIMATES OF 7decay 

Here, for the reader's convenience, we gather all the above estimates. 
Theorem 1 for m > 1 : 7decay > '^est', see (12) and (13); 

Theorem 1 for m = 1 : 7decay > (-Ist- 

Corollary 11 : 7dccay < r+ = min{(l + (/5fc)||6fc|| : ipk < 1}; 

7decay > = minA;(l - (/'A:)ll^fe|| if V^fc < 1 for ah k, 
where (pk = (^rj^)!^' h = B*Vk. 

Theorem 14 : 7decay ^ P foi' rn > 2, where p is defined in (36). 

We notice also that if these statements provide several lower or upper bounds for 7decay) then, 
obviously, one can take the best one of these. 



6. Numerical examples 

6.1. An example vi^ith 4 states and 2 controls. Take re = 4, m = 2 and consider the 
matrices 



A = i 



-b 








" 




' 1 


" 





—a 








B = 





1 




















a 





1 














b 







1 



where a, b are positive real numbers. Consider the LQR problem for (A, B) with Q = I , R = I . 
Table 1 collects the values of ||X|| and of 7dccay for different values of a, 6, obtained by numerical 
calculations. The last four columns of this Table show the values of r_, r_|_, £est and p, which 
are the lower and upper bounds for 7decay guaranteed by our theorems (see the previous section 
for a brief account). 

Row 1 shows that the bounds r_ , r+ for 7decay are very precise in the case of large separation 
of the spectrum of A. In rows 2 and 3, one can see that as the separation diminishes (and some 
ifkS approach to 1), the bounds r_, r+ become much more vague. 

In row 4, there is some k with ipk < 1, but we do not have (pk < 1 for all k. Hence, only the 
upper bound F-i- from Corollary 11 holds, and r_ is not defined. 

In rows 5, 6 and 7, > 1 for all k. Hence Corollary 11 provides no bounds at all, and we 
do not show the values of r_ , r_|_ . In these rows one can see how the lower estimate p for 7decay 
from Theorem 14 can give better results than iest from Theorem 1, especially if some eigenvalues 
of A are close together in comparison with 

Part (2) of Theorem 1 and Theorem 14 show that if the minimal singular value of B is large 
in comparison with the diameter of the spectrum of A, then the closed-loop spectrum divides 
in two parts: m eigenvalues are in the band Rez S ~"^crm] and the resting n — m 

eigenvalues lie in the band Rez e ( — \/3A, — max(£est, p)) • Within the values of a,b in the 
table, this result only applies to rows 6 and 7. For instance, for row 7, Part (2) of Theorem 1 
yields that two closed-loop eigenvalues lie in the band Rez £ [—1.4142, —0.8660] and two others 
in the band Rez G (-0.3811,-0.0198). Numerical simulation shows that two eigenvalues of 
A — BF satisfy Re 1^1,2 ~ —1.4024 and two others satisfy Re 1/3,4 w —0.1062. 

Simulation also shows that in many cases, the relative error in the estimate 7decay > ^Lt; 
which Theorem 1 gives for m = 1, is less than in the corresponding estimate for m > 1. (On 
the other hand, the quality of the control increases with the increase of m). 
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(X 


u 


II Y\\ 
W W 


Tdccay 


1 _ 


-L + 


f , 


P 


1 
i 


1 1^ 

iO 




1 ni 71 


n QQQQ 


u.yo ( u 


i.UioU 


n 70/1 n 


U.UoUD 


2 


5 


15 


1.0537 


9988 


0.8834 


1.1166 


0.6799 


0806 


3 


1.8 


5.4 


1.1667 


0.9910 


0.1006 


1.8994 


0.5403 


0.0806 


4 


1 


10 


1.1031 


0.9960 




1.1439 


0.3536 


0.0806 


5 


4 


4.1 


1.1456 


0.9928 






0.0018 


0.0806 


6 


0.2 


0.22 


5.4750 


0.2199 






7.07 • 10-5 


0.0396 


7 


0.1 


0.11 


10.2326 


0.1062 






1.77-10-5 


0.0198 



Table 1. Numerical results and bounds for the LQR problems, m = 2. 



6.2. A control problem for a mechanical problem. In many practical problems there is a 
large choice of possible physical or geometric configurations of the controller, which might make 
it necessary to solve a large amount of LQR optimization problems, in order to find a good one 
in some alternative sense. We will be speaking about the search of an LQR optimal regulator, 
which is also good in the sense that it has the largest possible 7decay 

In this subsection, we propose an algorithm which allows one to reduce drastically this search, 
by making use of our theoretical estimates. We will illustrate this algorithm on a simple me- 
chanical system. The same algorithm, in fact, can be applied to the following general class 
of problems: to optimize 7decay among a large finite family of LQR problems {A,Bj), with A 
skcw-Hermitian. In other words, the system matrix A is supposed to be fixed, but there are 
several possible choices for the control matrix B. 

This is not the only application of our bounds. We believe that in many cases the control 
designer can apply our results to obtain some a priori information on the systems in study. 

Consider a one-dimensional massless string. Attached to the string are N equal point masses 
of mass m, that are placed along it at equal distances h. It is assumed that the unperturbed 
string occupies the interval [0, {N + l)h] of the x axis in an xy plane; the string is supposed to 
move only in this plane. The two endpoints of the string are fixed, and it has constant tension 
T > 0. 

The problem is to stabilize the string using M controls, where 1 < M < TV. Namely, we 
choose point masses with numbers ji, j2; • • • ,jM, where 1 < Ji < J2 < • • • < Jm < -^j and 
apply a force to the point mass in the direction y. Every configuration (ji, j2, • • • ,Jm) of 
controls leads to its own linear quadratic control problem and to a corresponding stable closed- 
loop system, which is optimal in the linear quadratic sense. However, exponential decay rates of 
these closed-loop systems will depend on the chosen configurations of the control. The problem 
we discuss here is to find the configuration (ji, j2, ■ • • , Jm) which leads to the best exponential 
decay rate. 

In the experiment, we have chosen the parameters r/h = 10, m = 50 and N = 30. We tried 
the values M = 1,2,3,4,5,8. One can observe that 7dccay depends much on the choice of the 
configuration J' = {ji,j2, ■ ■ ■ iJm) (these are the numbers of the masses to which the control 
forces are applied). For example, if M = 5, then the best value of 7decay equals to 8.87 • 10"^, 
which is attained, for instance, for J = (2,5,11,19,27), while for J = (1,2,3,28,29,30) one 
only gets 7decay = 3.33 • 10^^, which is several times less. 

There are (^) configurations, and theoretically, the problem can be solved by a "brute force" 
complete search among all of them. However, even for moderate values of N and M, solving 
numerically (j^) LQR problems will be very time-consuming. 

If the position of the j-th point mass is {jh, yj), we obtain (in the linear approximation) the 
following system of ODEs: 

my'j = liVj+i + yj-i - 2%), j = 1,...,N, j ^ jk, 

' ^y'j = I + ~ 2%) +Uk, j =jk, k = l,...,M, 
^yo = VN+i = 0, 

where yo,yN+i have only been introduced for convenience in the notation. 
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where Cj S C is the j-th vector of the canonical basis. Then we obtain the control system 

x' = Ax + Bu. 

The energy of the system can be defined in terms of the following inner product in C^^: 
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The energy is E{x) = It is easy to show that energy is conserved, so that A is skew- 

Hermitian with respect to this inner product. 

Now we apply the Linear Quadratic Regulator using the cost functional 



l + \\u{t)fdt 



in order to stabilize the system. 

We can do a theoretical study of the system to obtain expressions to compute our estimates. 
Notice that our string is a very particular case of a nonhomogeneous string, whose spectral 
theory comes back to M.G. Krein, see [11, Section 8 of Chapter VI]. In our case, the eigenvalues 
of A are 

T17V 

sm I :rT^^ — 1 , -N <n<N,n^O, 



iXn 



-2i 



hm 



2{N + 1) 



and the corresponding orthonormal eigenvectors are where 

nkir 



sm 



N + 1 



l<k<N 



See the paper [23] by Micu, where the same matrix A appeared in the context of a semidiscrete 
numerical scheme for ID wave equation. We also refer to [1], [24] for a related inverse problem. 

The operator \/N + 1 Bq maps the canonical basis of C*^ onto an orthonormal system of M 
vectors in C'^^ (we use the inner product (•, ■)e in and the standard one in C^). Hence, 
is an isometry and it follows that 

\\B\\=aM{B) = —^. 

m\/ N + 1 

Finally, the vectors bk = B*Vk can be computed to obtain 



M 



m2(Ar + 1)2 ^■'''^ + 1 



sm 



-N < k < N,k y^O. 



Using Corollary 11, we can give an upper bound for 7decayi assuming that some (pk > 1- 
Theorem 1 and Corollary 1 1 (if it applies) can be used to obtain a lower bound for 7decay • The 
following algorithm uses these bounds to reduce the number of LQR problems being computed. 
In the course of its execution, the upper and the lower theoretical bounds for all configurations 
are taken into account, but the LQ optimal regulator is actually computed for a fewer number 
of configurations. 

The algorithm works as follows: 

(1) Calculate the eigenvalues iXk and the corresponding eigenvectors Vk of A. 
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M 


Tdecay 


Time (s) 


LQRs computed 


% computed 


1 


6.53 • 10"^ 


2 


30 


100 


2 


3.39 • 10"^ 


4 


76 


17.5 


3 


5.82 • 10-"^ 


28 


441 


10.9 


4 


7.70 • 10-^ 


123 


481 


1.76 


5 


8.87-10-^ 


738 


5505 


3.86 


8 


1.2 • 10-3 


56675 


198369 


3.39 



Table 2. Results for N = 30,T/h = 10, m = 50. 



(2) For each control configuration J = (ii, . . . , Jm), compute the vectors bk and the quan- 
tities Uj and Lj, which are the upper and the lower theoretical bounds for 7decay Set 
Uj = +00 if an upper bound is not available. 

(3) Select the configuration having the maximal L j. Solve the LQR problem numerically 
for this configuration and compute 7decay 

(4) Now we proceed to a search, defined recursively as follows. Let 7 be the best 7decay found 
so far. If for all configurations J whose corresponding 7decay has not been computed 
yet, Uj < "y, the search stops, and this current value of 7 is taken for the optimal 7decay 
If there are configurations whose 7decay has not been computed that have Uj > 7, 
the algorithm selects the one having the greatest Uj. For this configuration, it solves 
the LQR problem numerically, computes its 7decay and updates 7 according to the rule 
7 := max(7, 7decay)- This is the best 7decay found so far. 

(5) The algorithm stops after having exhausted all possible configurations. It returns the 
last value of 7, which is equal to the maximum of the values of 7decay over all possible 
configurations. 

Observe that this algorithm also allows one to compute all the configurations having the 

optimal 7dccay 

The results of the execution of the algorithm are shown on Table 2. The computations were 
done on a modern desktop computer. Recall that we have chosen the total number of masses 

= 30. The table shows that the decay rate 7decay improves when M increases. The fourth 
column collects the number of LQRs the algorithm had to solve, and the fifth column shows the 
ratio between the total of (^) possible configurations and the number of configurations that 
were actually processed. One can see that in many cases, our algorithm reduces drastically the 
amount of computations. 

The values r/h = 10 and m = 50 have been chosen for these computations because they 
provide a moderate separation of the spectrum of A with respect to If we fix m = 50 

and increase r/h (say r/h = 1000), then the number of computations is further reduced, since 
the separation of the spectrum of A increases and we obtain tighter theoretical bounds. On the 
other hand, if one sets r/h to a small enough value while maintaining m fixed, our algorithm 
will not provide much save in the computations. 



7. Some open questions 

Question 1. Assume that m<n, R = I,Q = I and that a skew-Hermitian matrix A is 
fixed. Does it follow that there is a constant C = C{A) such that 7dccay < C, independently of 
B? Part (2) of Theorem 1 shows that it is so if, for instance, am{B) > or even if we assume 

that (Jm{B) > where / is any function on [0, +00) such that lim^^-^oo f{x) = +00. We 

conjecture that it is true in general. 

Question 2. Corollary 8 can be used to obtain an upper bound for 7decay However, if (pk < 1 
for some k, then Corollary 11 gives a much tighter bound. Can one give a tighter upper bound 
even when ipk > ^ for every k? 
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8. Conclusions 

• The bounds iest-, (-\st given in Theorem 1 can be applied only if all the eigenvalues of A 
are different. 

• The lower bound p given in Theorem 14 is the one which can be used in a more general 
setting (namely > 0, which allows some eigenvalues of A to coincide). There are 
cases when it is the best bound available. It happens, in particular, if some eigenvalues 
of A are close together (compared with ||-B||). 

• The two-sided bound given in Corollary 11 holds only when (pk < 1 for all A;, i.e., when 
the spectrum of A is separated enough. 

• If all (pk are small, this two-sided bound is very tight and one can take do = min ||5fc|| as 
a good approximation for 7decay 

• When all are small, one can also use Theorem 10 to locate with precision all the 
eigenvalues of the closed-loop system. 

• Corollary 8 shows that if m < n and the diameter A of the spectrum of A is much 
smaller than all singular values of B, then 7decay is less than -y/S A. 

• One can observe that, as a rule, if the separation of the eigenvalues of A increases or the 
number of controls m increases, then 7decay grows. 

• If one has to find an optimal 7decay among a large finite family of LQR control problems, 
our estimates permit one to design an algorithm to reduce the search (in some situations, 
drastically; see Subsection 6.2). 

• By now, we only have estimates of 7decay for the case of a skew-Hermitian matrix A. It 
would be very desirable to generalize our main results to non-skew Hermitian matrices, 
or at least to the case of matrices A such that Re A > 0. 
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